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Abstract: The Gerda experiment searches for the neutrinoless double beta decay of ^^Ge us- 
ing high-purity germanium detectors enriched in ^^Ge. The analysis of the signal time structure 
provides a powerful tool to identify neutrinoless double beta decay events and to discriminate 
them from gamma-ray induced backgrounds. Enhanced pulse shape discrimination capabilities of 
Broad Energy Germanium detectors with a small read-out electrode have been recently reported. 
This paper describes the full simulation of the response of such a detector, including the Monte 
Carlo modeling of radiation interaction and subsequent signal shape calculation. A pulse shape 
discrimination method based on the ratio between the maximum current signal amplitude and the 
event energy applied to the simulated data shows quantitative agreement with the experimental data 
acquired with calibration sources. The simulation has been used to study the survival probabilities 
of the decays which occur inside the detector volume and are difficult to assess experimentally. 
Such internal decay events are produced by the cosmogenic radio-isotopes ^^Ge and ^*^Co and the 
neutrinoless double beta decay of ^^Ge. Fixing the experimental acceptance of the double escape 
peak of the 2.614 MeV photon to 90%, the estimated survival probabilities at Q^^ = 2.039 MeV 
are (86 ± 3)% for ''^Ge neutrinoless double beta decays, (4.5 ± 0.3)% for the *'^Ge daughter *'^Ga, 
and (0.9^°;^)% for ^^Co decays. 
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1. Introduction 

The search for neutrinoless double beta decay (Ovj8j8) is the only practical method to probe whether 
the neutrino is a Majorana particle, i.e. its own anti-paiticle. Ovj3/3 decay violates lepton num- 
ber conservation and would establish new physics beyond the Standard Model of particle physics. 
Provided that the exchange of light Majorana neutrinos are generating the leading term in OvjS j3 de- 
cay, effective neutrino masses down to a few tens of meV will be probed in the next generation of 
Ovj3j3 decay experiments. 

One promising isotope is ^^Ge. High-Purity Germanium (HPGe) detectors can be produced 
from germanium material enriched in '^Ge, which serves simultaneously as a source and as detector 
medium. The advantages of the intrinsic radio-purity and excellent spectroscopic performance 
(2 — 3 keV full width at half maximum resolution Qpp = 2039 keV) of HPGe detectors have 
been recognized early [HU, and until today, the most stringent experimental limits on Ovj3j3 decay 
have been achieved with germanium crystals enriched in ^^Ge [0]. Also one group reported a claim 
of evidence [2D- 
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The research work presented here was carried out in the framework of the GERmanium 
Detector Array (Gerda) experiment O] which recently started its commissioning phase at the 
Laboratori Nazionali del Gran Sasso (LNGS) of the INFN in Italy. An array of bare HPGe 
detectors enriched to 86% in ^^Ge will be immersed in a cryogenic liquid which serves simul- 
taneously as a coolant for the germanium detectors and as a high-purity shield against exter- 
nal radiation. The experiment pursues a staged implementation: the mentioned claim of evi- 
dence will be scrutinized in its first phase with about 15 kg • years exposure and a background 
index at Qp^ of < 10^^ counts/(keV • kg • year). The second phase of the experiment will ex- 
plore half-lives up to 2 • 10^^ years, with 100 kg • years of exposure and background index of 
< 10^-^ counts/(keV • kg • year). Contingent on the results of the first two phases, a third phase 
is conceived to probe half-lives > 10^^ years. The corresponding effective neutrino mass of 
> 10 meV is predicted based on the results of neutrino oscillation experiments assuming an in- 
verted mass hierarchy To explore this parameter regime, an exposure of several 1000 kg • years 
and backgrounds < 10^^ counts/(keV • kg • year) ai^e required. To reach such background levels, 
which are about three orders of magnitude below the best current values, novel techniques ai^e 
required which exploit the decay characteristics and topology of background events. 

In recent works, we could show that p-type HPGe detectors produced by Canberra Semicon- 
ductor 10] with a small read-out electrode, referred to as thick-window Broad Energy Germa- 
nium (BEGe) detectors, exhibit pulse shape discrimination performance [0, superior to coaxial 
HPGe detectors. This allows to efficiently distinguish between single-site events (OVj8j8-like) 
and multi-site events (gamma-ray backgrounds). Investigations of a bare BEGe detector in liquid 
argon showed unaltered performance as compared to operation in a vacuum cryostat [0]. An ad- 
ditional benefit of the small size of the read-out electrode is its low capacitance which results in 
a lower noise and therefore a superior resolution and a lower threshold compared to other types 
of large-volume HPGe detectors. Following a full production chain test f lTOl of BEGe detectors 
with modified isotopic composition (depleted in ^^Ge), the Gerda collaboration decided to adopt 
thick-window BEGe detectors as their baseline design for the second phase of Gerda. R&D with 
BEGe detectors is also carried out by the Majorana collaboration flTTl IT^ . 

The scope of this paper is the ab-initio modeling of the signal pulse shapes of BEGe detectors 
created by particle interactions and the comparison with experimental data. The energy depositions 
and event topology were generated with the Geant4 simulation package [1T3L IHH and the conse- 
quent signal pulse shape generated by drifting the charge carrier in the calculated electric field. The 
characteristic pulse shapes of the simulated pulses are then analyzed. First we validate the simula- 
tion of pulse shapes and discrimination methods with experimental data using calibration sources 
located externally to the crystal. Then, we calculate the expected pulse shape discrimination (PSD) 
cut efficiencies for ^°Co, ^^Ga and OvjSjS decays internal to the detector. We also study the pulse 
shapes of energy depositions close to the small read-out electrode and in its vicinity. 

2. Overview of the simulation 

The simulation presented here replicates the physical processes involved in the generation of the 
Ge detector signals with application to the particular case of Broad Energy Germanium (BEGe) 
detectors. The simulation can be divided into three logical blocks: 1) interaction of ionizing par- 
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Figure 1. Data flow of the simulation organized into three logical blocks. The data flow includes a first phase 
(thin arrows) in which a library of pulses and baselines is created and a second phase (thick arrows) in which 
the libraries are used to generate the detector signals. In block III, FE & DAQ refers to the preamplifier and 
the sampling device. 

tides with the Ge crystal; 2) collection of the charge carriers produced at each interaction point 
at the detector electrodes and time evolution of the induced signals; 3) signal shaping by read-out 
electronic devices. 

The first block consists of a Monte Carlo simulation, including all the physical processes 
involved in the passage of gamma-rays or charged particles through matter. It provides the inter- 
action points and the corresponding energy losses within the crystal. This part of the simulation 
is performed within the MaGe framework [1T3I] . which is based on the Geant4 simulation pack- 
age [D, 11. 

The second block describes the dynamics of the charge carriers generated inside the detector 
and provides the signals induced on the electrodes by the charge movement. It is calculated by 
using a modified version of the Multi Geometry Simulation (MGS) software flTOl . 

The last block of the simulation takes into account the effects of the read-out electronic devices 
and of the electronic noise. It generates signals that can be directly compared to the measured ones. 

Figure |T] illustrates the implementation diagram of these logical blocks in the software. In the 
following two subsections the second and the third block of the simulation are explained in more 
detail. 

2.1 Field calculation, charge carrier transport and signal generation 

The pulse simulation performed with the MGS-based software includes the following steps: 1) com- 
putation of the electric field inside the detector at defined potentials of the electrodes; 2) transport 
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of charge carries toward the electrodes; 3) time evolution of the signal induced at the read-out 
electrode by the moving charges. 

Since a semiconductor detector can be considered as an electrostatic system, the electric field 
can be computed by solving the Maxwell equations or, equivalently, by solving the Poisson equa- 
tion with a set of boundary conditions for the potential. In p-type Ge detectors, the p-n junction 
forms at the border of the ~ 1 mm thick donor-doped surface layer (n-i- layer). At operational bias 
voltages the p-type volume is fully depleted of free charges (active volume), while most of the n+ 
layer retains electrons in the conduction band (anode, forming a dead layer). In our simulation we 
approximate this situation by defining the active volume with a negative space charge distribution 
p proportional to the net p-type impurity distribution, adjoined to the dead layer on which the bias 
voltage is applied ^ The boundary conditions are then provided by the value of the potential on the 
conductive surfaces: Vcathode and Vanode- Therefore the system of equations to solve is: 

' V^<^{r) = -p{r)/{eoer) 

' ^\Scaihode = ^cathode (2-1) 
'i'\s anode ~ ^aiiode 

where r is the position, if) is the electric potential, p (r) is the charge density distribution, £o is 
the vacuum electrical permittivity, ~ 16 is the relative permittivity of Ge, <^\s is the potential at 
the surface S surrounding the considered electrode. 

The movement of charge carriers within the active volume is computed by using two phe- 
nomenological models lUTl ITEH which provide the drift velocity as a function of the electric field 
magnitude and of the electric field direction relative to the crystallographic axes. We assume that 
each transfer of energy to the Ge crystal lattice results in the generation of a cloud of free electrons 
and holes, which further drift as two independent clusters. These ai^e approximated in our simu- 
lation by two point-like charges with opposite sign. The trajectories are calculated with a fourth 
order Runge-Kutta method with 1 ns time step. 

Using the simulated trajectories of charge carriers and the weighting potential^ distribution 
inside the detector, the charge signal Q{f{t)) induced on the electrodes is computed by using the 
Shockley-Ramo theorem []T3 : 

Q{f{t)) = -qtot^.ir{t)) (2-2) 

where r{t) is the position of the charge cluster at the time t, qtot the total charge of the cluster and 
^{r{t)) is the weighting potential. 

The total signal of a simulated particle event composed of several interactions is calculated as 
a sum of pulses from each hit, whose position and energy deposition are provided beforehand with 
the MaGe Monte Carlo simulation. To reduce the computation time, the generation of the signals 
corresponding to the individual hits is performed by using a library of pre-calculated pulses. The 
library is generated by dividing the detector active volume in 1 mm^ cubic cells, and simulating 

'since our field calculation is performed on a 0.5 mm step grid, it can not describe the thin part of the n+ volume in 
which the conduction-band electrons are depleted (the n-side of the p-n junction, with positive space charge). 

^The weighting potential is a dimension-less quantity defined as the electric potential calculated when the considered 
electrode is kept at a unit potential, all other electrodes are grounded and all charges inside the device are removed. 
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a single normalized interaction at each comer. The library must be generated only once for each 
simulated detector geometry and bias voltage setting. Then each Monte Carlo generated interaction 
is associated to one of the cubic cells in the pulse library. The signal is computed as the weighted 
average of the eight pulses in each comer of the cube, where the weight is given by the inverse of the 
cubic euclidean distance between the interaction point and the considered corner. The amplitudes 
of the individual interaction pulses are then normalized according to the energy depositions in the 
hits and all the pulses of the event are added up to one combined signal. 

The use of a pulse library decreases significantly the processing time of the simulation in such 
a way that the pulse computation time is comparable with that of the Monte Carlo simulation. 

2.2 Read-out electronics response and noise 

The read-out electronics includes the preamplifier and the digital sampling device. The response 
function was determined by providing the preamplifier input with an impulse generated by a high- 
precision pulser, and then by deconvolving the digitally sampled signal with the input signal. The 
simulated detector signals are consequently convolved with this response function. This is per- 
formed for each pulse in the library mentioned in previous subsection, so that it needs to be done 
only once for each simulated experimental setup. 

To reproduce the electronic noise present in the experimental data, samples are taken from a 
library of experimentally recorded baselines. The amplitude of the noise is normalized according 
to the experimental signal-to-noise ratio, and the noise sample is then added to the calculated full 
event signal. The output of the simulation is a file of signals similar to those recorded experimen- 
tally with a digital data acquisition system. Therefore it is possible to apply the same analysis tools 
to both experimental and simulated data. 

3. Modeling of BEGe detectors 

The detector used for the present simulation is a modified thick-window version of a Broad Energy 
Germanium (BEGe) detector [EOl . a standard p-type High-Purity Ge (HPGe) detector offered by 
Canberra Semiconductor. These detectors have a cylindrical shape and a small B-implanted p-i- 
electrode on one of the flat surfaces. The Li-diffused n-i- electrode (between 0.4 mm and 0.8 mm 
in the "thick window" modification) covers the rest of the outer surface and is separated from the 
p-H electrode by a circular groove. The detector is mounted in a 1.5 mm thick cylindrical aluminum 
housing. The positive high voltage is applied to the n-i- electrode while the signal is read out from 
the p-i- electrode where the holes are collected. For the purpose of the present work the specific 
case of the BE3830/S model was considered. The BE3830/S model has a diameter of 71 mm and a 
thickness of 32 mm. The detector was experimentally characterized in Refs. lETl E3l ■ A schematic 
drawing of a typical BEGe detector is shown in Figure 0. 

In this section the signal formation and development in the BEGe detectors will be discussed 
following the steps outlined in Section PTT] . Moreover, the results of the simulation will be used to 
show how it is possible to distinguish between single-site and multiple-site events by analyzing 
the pulse shapes. 
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n+ electrode 
(high voltage contact) 



p+ electrode (signal read-out contact) 



Figure 2. Sketch of a BEGe detector. The signal read-out electrode and the groove are not to scale. The 
plot shows a vertical section of the detector passing through the symmetry axis. 



3.1 The electric field inside BEGe detectors 

The pulse shape discrimination properties of BEGe detectors have been attributed to the peculiar 
internal electric field which is created by the small size of the read-out electrode. In Section ^71] 
the space charge distribution inside the active volume (equivalent to the net distribution of acceptor 
impurities) was identified along with the electrode potentials as a source of electric field inside the 
detector. To better understand the two contributions, the linear superposition principle can be used 
to separate the potential into the two individual components: 



<^(r) = (/)o(r) + (/)p(r) 



(3.1) 



where <po is the potential calculated considering only the electrode potentials and no impurity charge 
(p (r) = Vr) and <pp is the potential generated by the impurity charge distribution when ground- 
ing all the electrodes. Therefore, we can solve the following two problems and then add up the 
solutions: 

I Scathode 
0pl 



*/*0 I Scathode 
S anode 





^cathode 
^anode 



-p{f)/ieoer) 








(3.2) 



Similarly, since the electric field is determined by the linear relation E 
separated into two components: 

E{r)=Eo{r)+Ep{r) 



-V(j), it can be also 



(3.3) 



where Eo(r) = -V^oir) and Ep{r) = -V0p(r). 

Figure shows the electric potential and the electric field strength of the two contributions and 
their sum for the BE3830/s detector operated in its nominal configuration, i.e. cathode grounded, 
anode set at 3500 V and p ~ 10^° impurity atoms/cm^. 

Figure 0.a and ^b show the electric potential and field strength generated only by the elec- 
trodes (00 (r) and | |£o C?^) 1 1 )■ The potential and the field show a sharp variation in the region close to 
the small-size p-i- electrode. The field in the rest of the volume is so weak that the collection time 
would be longer than the characteristic recombination time and most of the charge would be lost. 

Figure 0.c and 0.d show the electric potential and field strength provided only by the impurity 
charge distribution (0p(r) and ||£p(r)||). For the considered BEGe detector the charge concentra- 
tion is approximated by a uniform distribution. The resulting negative electric potential reaches its 



-6- 




(a) M7) 



4.0e3 



3.()e3 



2.0e3 -3 



l.()e3 a. 



O.OeO 



1.2e6 




O.OeO 




O.OeO 



-5.0e2 



-1.0e3 -S 



-1.5e3 a. 



-2.0e3 



1 



4.0e3 



3.0e3 



(c) (ppir) 



(e) (j,{T-) 



3.0e5 




2.0e3 -3 



1.0e3 a. 



O.OeO 



1.6e6 




(b) \\Eo{r) 



(d) l|£p(?) 



Figure 3. Simulated electric potential and electric field strength for different configurations of a BEGe 
detector In (a) and (b) the electrode potential is considered, in (c) and (d) the charge distribution, and in (e) 
and (f) the sum of the two contributions. The plots show half of a vertical section of the detector passing 
through the symmetry axis. The cathode is drawn in red and the anode in black. 



peak value in the middle of the detector. The direction of the electric field in the region close to the 
small cathode is opposite to that of Eq, whereas in the rest of the detector the field helps to move 
the charges produced close to the outer n+ electrode towards the central slice of the detector. 

Figure ^e and ^.f show the total electric potential and field strength (0 (r) and | \E{r) 1 1). The 
potential close to the small p+ electrode is dominated by the electrodes contribution smoothed out 
by the opposite contribution of the impurity charge distribution field. In the rest of the volume the 
dominant contribution is provided by Bp. The effect of the Ep field is to bring the holes in the 
center of the detector while the field Eq subsequently collects them to the read-out electrode. This 
peculiar way of charge transportation in BEGe detectors leads to the favorable signal shape as it 
will be discussed in the next sections. 

3.2 Signal development in dependence of the interaction position 

As given by equation (|2.2|) , the weighting potential determines the charge that a cluster of charge 
carriers drifting inside the detector induces on the read-out electrode. Figure ^ shows the weighting 
potential (j)w{r) and the strength of the weighting field ||£vy(r)|| = || — V0,i.(r)|| of the cathode in 
a BEGe detector. The weighting potential (similarly to 0o in Figure ^a) has a sharp variation 
in the region close to the small-size p+ electrode and it is very weak in the rest of the detector 
volume (blue area in Figure If an interaction occurs in the volume of weak (p^, according to 
equation (|Z.2| ), the signal induced by the holes drifting towards the small cathode remains small 
until the charges arrive at about 1 cm away from the electrode and then rapidly grows until the 
holes are collected. 
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(a) K(r) 

Figure 4. Weighting potential (a) and weighting field 
computed for half of a vertical section of the BEGe 
(density maps) and for the symmetry axis (black lines 
anode in black. 




strength (b) of the small read-out electrode (cathode) 
detector passing through the detector symmetry axis 
in the side plots). The cathode is drawn in red and the 



On the other hand, in most of the active volume the electron contribution will be present only 
at the beginning of the pulse, since electrons are collected to the outer n+ electrode and the portion 
of the detector volume that is close to the outer surface is much larger than the portion close to 
the p+ electrode. This effect is further enhanced due to the fact that electrons move with velocities 
roughly two times higher than the holes. Moreover, the signal induced by the electron drifting will 
be relatively small, as most electron clusters move in the region of weak . Thus, the contribution 
of the electrons to the formation of the signals is expected to be negligible for most of the interaction 
positions. 

Figure ^.a shows the electron and hole trajectories for three interactions in the "bulk" detector 
volume far from the p+ electrode. The holes follow the electric field (Figure ^f) and are first 
transported into the middle slice of the detector, then drift towards the center of the detector and 
finally their trajectory bends towards the read-out electrode. It can be noticed that for all these 
events, the last part of the hole collection happens along a common path which is independent of 
the starting position. Accordingly, the last part of the induced signal is identical for the different 
events. Since the first part of the signal, induced by the holes and electrons in weak regions, is 
comparatively small, the signal shapes are essentially independent of the interaction position. This 
can be clearly seen in Figure ^b, which shows the charge and current pulses corresponding to the 
trajectories from Figure ^a. The only visible difference is the time shift of the rising part, which 
depends on the length of the charge carrier path to reach the strong (p^ region. It is important to 
stress that this type of almost indistinguishable signals originates from interactions in most of the 
detector volume, including corner regions. For further discussion we will refer to such signals as 
"type I" trajectories. 

The shape of the current signals in Figure ^.b can be understood by differentiating the function 
in equation (|2.2|) . The current induced at the cathode by a charge carrier is then given by: 

nm) = = qtot m))-E^m) (3.4) 
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Figure 5. Simulated electron-hole trajectories (a) and corresponding charge and current pulses (b) for 
three events occurring in different places in the bulk detector volume far from the p+ electrode ("type I" 
trajectories). The small oscillation in the current signals after the peak originates from the experimentally 
measured FE and DAQ response the pulses are convolved with. 



From this equation we can see that the induced current I{f{t)) depends on the velocity v{f{t)) and 
the weighting field Ew{r{t)) at the position f{t) of the charge cluster. The charge carrier velocity 
can vary roughly between 5 cm/s and 10 cm/s [ETH . while E„. can increase by more than a factor 
of 20 close to the read-out electrode (see Figure ^.b). It is evident that E„ has a dominant effect 
on the current signal - the signals of "type I" trajectories in Figure ^b feature a significant current 
peak at the end of the hole collection. 

Two other types of less common trajectories can be identified when the interactions occur in 
the close vicinity of the groove and the small p+ electrode. Examples of these kind of events are 
displayed in Figure ^ The "type 11" trajectories originate close to the p+ electrode (green and black 
color). In these events the holes are directly and quickly collected at the cathode. Now also the 
electrons drifting in the opposite direction provide a significant contribution to the signal, since they 
aie now moving in a region of strong (p^^.. The closer the interaction occurs to the cathode, the more 
important the signal induced by electrons becomes. The signal is fully dominated by the electron 
contribution for interactions within ~ 2 mm from the p+ electrode (black example in Figure The 
induced charge signal rises quickly at the beginning and then, as they drift away from the cathode 
into the weaker <p„ regions, the signal growth slows down. The current peak appears at the very 
beginning of the collection time. For events occurring few mm to ~ 1 cm from the p+ electrode 
(green example in Figure ^ neither electrons nor holes traverse the full thickness of the strong ^„ 
region. So the main part of the signal is induced in a relatively short time and the rise time is thus 
faster. The current peak is amplified because contributions from both charge carrier types add up. 
The current amplification can be further enhanced if the interaction happens close to the inner edge 
of the groove, because here is strongest (Figure ^.b). 

Starting points in a zone close to the anode, ~ 1 .5 cm from the detector symmetry axis, result in 
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(a) Trajectories (b) Charge and current pulses 

Figure 6. Simulated electron-hole trajectories (a) and corresponding charge and current pulses (b) for three 
events occurring next to the cathode and a "type I" event (orange line). For the green and black events ("type 
11" trajectories) both the electron and hole contribution are present at the very beginning of the signal. In the 
blue event ("type III") the electron collection is quick and provides a characteristic kink in the first part of 
the signal. 




"type III" trajectories (shown blue in Figure p. For these events the electrons are collected quickly, 
and since (p^ is still noticeable in this region, they provide a characteristic kink in the first part of 
the signal. This quick increase at the beginning causes the 10% to 90% rise time measurement to 
give higher values for these events than for the "standard" type I trajectory events. 

The three types of signals and the extent of the volumes from which they originate can be 
better understood from the plots shown in Figure [7[ The figure shows the rise times in the 10% to 
30% and 10% to 90% intervals of the signal height, as a function of the interaction position. The 
type II trajectories originate in the volume close to the p+ electrode (~6% of the active volume), 
distinguished by short rise times plus a small zone of longer rise times few mm away from the center 
of the electrode. The type III trajectory starting points can be identified in the area of increased rise 
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Figure 8. Distribution of the maximum current pulse amplitude A for simulated single interactions with unit 
energy as a function of the interaction point. A is constant in most of the detector volume (type I trajectories), 
but it is amplified in the region close to the cathode (type 11 trajectories). The region of lower A values close 
to the outer radius of the groove is an artifact due to the use of the pulse shape library with 1 mm step size. 
Here the charge collection time varies on scales smaller than the library grid with the effect that the averaging 
of library pulses leads to a reduction in current-peak amplitude. The effect is not present when the signals 
are generated directly without the use of the hbrary. 

time beyond the outer diameter of the groove (~1% of the total active volume). Interactions in the 
rest of the detector volume result in the most common type I trajectories. 

3.3 Discrimination between single-site and multiple-site events 

In the previous subsection we have shown that most of the single interactions in BEGe detectors 
produce uniform signals (type I trajectories), differing one from the other only by the total charge 
collection times. In addition, the current signals have a simple shape with only a single narrow peak 
at the end of charge collection. These features can be exploited for a powerful discrimination of 
Single-Site Events (SSE) and Multiple-Site Events (MSE). A discrimination method based on the 
current signal amplitude was introduced in Ref. |]7|] and explained by using an empirical estimate of 
the weighting fields in BEGe detectors flSD. Here we recall the basic idea of the method and refine 
the discussion based on more accurate electric potential and field calculations. 

For typical BEGe events the part of the hole trajectory that passes through strong Eyt- is al- 
ways the same. Therefore, according to equation (|3.4|) the amplitude of the induced current signal 
depends only on the total charge of the considered hole cluster. Figure |S| shows the value of the 
maximum current pulse amplitude A for simulated single interactions with unit energy deposition, 
as a function of the interaction position in a vertical section of a BEGe detector. We see that the 
parameter A is constant in most of the detector volume (corresponding to the typical type I trajec- 
tories) except for the region close to the p-i- electrode, where the type II trajectories with amplified 
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Figure 9. The upper panel shows typical charge pulses for SSE and MSE while the lower panel shows the 
corresponding current pulses obtained as the derivatives of the charge signals. The dashed lines in the top 
panel show the contributions of the individual single interactions to the total charge pulse in the MSE. 



current signals originate (Figure p. 

The Pulse Shape Discrimination (PSD) method uses the parameter A normalized to the to- 
tal event energy E: the AjE ratio. The concept is depicted in Figure ^ For SSE, all energy is 
transferred to a single charge cluster - it follows that AjE \& constant for such events. In MSE, 
the total event energy is shared between several spatially-separated charge clusters (Figure ^ illus- 
trates the case of three clusters). Since the cuiTcnt peaks are naiTow and the charge collection time 
is position-dependent, the AjE ratio will be smaller than the constant value given by SSE. 

The region of amplified cuiTcnt signals (approximately a hemisphere of a ~ 13 mm radius 
from the center of the p-i- electrode in Figure ^ coiTcsponding to ~ 5% of the detector active 
volume) was already identified in Refs. [0, ED- It limits the efficiency of our PSD method, since 
interactions from MSE occurring in this volume can have their A /£ ratio amplified above the SSE 
discrimination threshold. On the positive side, this effect can be used to identify surface events 
occurring on the p-H electrode and groove surfaces. 

4. Validation of the simulation 

The implementation of BEGe detectors in the MaGe Monte Carlo simulation framework, accu- 
rately reproducing their radiation detection efficiency, was presented already in [EL EL El, E2l- Here 
we report only the measurements carried out to validate the pulse shape simulation. Two sets of 
measurements were performed for this purpose. First we used a collimated ^^^Am source to gen- 
erate well localized interactions and to compare directly the pulse shapes for different interaction 
positions close to the surface of the detector. Then, a ^^^Th source was used to investigate the 
distribution of the pulse rise times and of the parameter A/E" as a function of the energy for events 
in the whole detector volume. 
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4.1 The experimental setup 

The detector used for the vahdation measurements was the BE3830/s BEGe detector described in 
Section 0. The front end read out of the signals was performed with the Canberra charge sensitive 
preamplifier 2002CSL [EOl . which is integrated in the housing of the detector. The preamplifier 
output was digitally recorded with a 4 channel N1728B CAEN NIM flash analogue-digital con- 
verter [E3I1 running at 100 MHz sampling frequency with a precision of 14 bits. This module is 
fitted with a USB connection for communication with a PC. To control the digitizer setup, the 
acquisition parameters and the storage of the data, the PC was interfaced to the NIM module by 
using the TUC software \UM . We recorded pulse shapes with a total length of 40 /is including a 
baseline before the signal of ~ 10 /is. The energy reconstruction was performed off-line by using 
the GELATIO software [E5l1 implementing the Moving Window Deconvolution approach lESII and 
several digital filters were applied to the data to reject pile-up and noise events. 

To perform the ^^'Am measurements, a mechanical device was built to allow the movement 
of coUimated sources with sufficient accuracy along the diameter of the front face and circularly 
ai^ound the symmetry axis of the detector. The collimator had a hole of 1 mm diameter and a length 
of 34 mm. These ^^^Am coUimated measurements were accurately simulated by using a conical 
beam of 59.5 keV photons. 

4.2 Pulse shape comparison with low energy gamma-ray beams 

Using a coUimated beam of low energy gamma-rays (59.5 keV) from a Am source allows to ob- 
tain pulse shape data with well defined interaction coordinates. As the 59.5 keV photons penetrate 
only a few millimeters underneath the germanium crystal surface and the volume spread of their 
energy deposition is of similar size, the narrow beam can be used to study pulse shape position de- 
pendence with few millimeters resolution. As the topologies of the individual photon interactions 
are statistically variable, we selected events with energy deposition corresponding to the 59.5 keV 
peak and calculated the average signals for each position of the coUimated source. This averaging 
procedure also reduced significantly the electronic noise. 

The pulse shape comparison included scanning along the diameter of the detector and a cir- 
cular scan at a fixed distance from the symmetry axis of the detector. The radial scan allows the 
study of pulse shapes as a function of the distance from the center. With the circular scan the 
rise time variations due to the drift velocity anisotropy caused by the crystal lattice structure were 
investigated. 

The radial scan was performed by moving the coUimated ^^'Am source along the radius of 
the detector on the top surface of the end-cap in steps of 0.5 cm. The simulated hole trajectories 
for each collimator position and the average signals are shown in Figure [T^. Good agreement is 
observed between the simulated and experimental pulses. 

As an increase of the radius value corresponds to an increment of the hole collection time, 
the total rise time should increase as the interaction point moves away from the symmetry axis. 
However, according to the discussion in Section for interactions far from the p-i- electrode the 
pulse shapes are expected to be different only at the very beginning. Looking at Figure [7[b BEGe 
signals along the top surface show only minimal variation of the 10% — 90% rise time. In order to 
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Figure 10. Simulated hole trajectories (a) and average signals (b) computed for several interaction points 
at different radii on the detector top surface. The interactions were obtained by moving a 1 mm collimated 
^'^^ Am source on the top surface of the end-cap along the radial direction. The small deviation present at the 
very beginning of the signals can be related to the fixed grid width used for defining the library of pulses. 
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Figure 11. 1% — 90% and 10% — 90% rise time values reconstructed from the simulated and experimental 
average pulses for different interaction points. The interactions were obtained by moving a 1 mm collimated 
^^'Am source on the top surface of the end-cap, along the radial direction (a) and around a part of a circle 
centered on the detector symmetry axis (b). The rise time uncertainty was estimated by determining the 
spread of rise time values computed from several measurements at one source position. 



observe a significant effect, the rise time between 1% and 90% has to be considered. This can be 
clearly seen in Figure [TTj-a for both the experimental and simulated data. 

For the circular scan the collimated Am source was moved around a circle with a radius of 
34 mm centered on the detector's symmetry axis on the top surface of the end-cap with a ~ 10° 



-14- 



step. The drift velocity anisotropy arises due to the charge carriers drifting at different angles 
relative to the detector's crystallographic axes. Like with the radial scanning, also in this case the 
differences occur only in the first part of the charge collection, and the 1 % — 90% rise time has to 
be used. The comparison of circular-scan rise times as a function of the circumference angle from 
the simulation and the experiment is shown in Figure PT|.b. 

Although the experimental data show a behavior coherent with the simulation in both the radial 
and circular scans, the agreement is only qualitative. The 1% — 90% rise times of the simulated 
data are systematically shifted by 20 — 40 ns (5% — 10% of the total rise time) to lower values than 
the experimental data. In contrast, the 10% — 90% rise times of the simulated data are shifted by 
20 ns to higher values than the experimental data. The change in the sign of the discrepancy is due to 
differences in the shape of the average pulses. Furthermore, the oscillation amplitude due to crystal 
anisotropy in Figure [rT].b is higher in the simulation than in the experimental data by 30%. Unlike 
the experimental data, the simulation shows a residual oscillation of about 5 ns also in the 10% — 
90% rise time plot. From these comparisons we can conclude that the simulation is more sensitive 
to the drift velocity anisotropy than the real detector and probably further small discrepancies in 
the signal calculation are present. This allows us to identify several areas of possible improvements 
of our simulation, including the accuracy in the geometry definition (especially important in case 
of the p+ electrode), the finite grid width used for defining the library of pulses and the charge 
carrier mobility model parameterizations (defined by using a different crystal [ITElD . Also further 
measurements could be performed to characterize the position dependence of the BEGe signals 
(e.g. by performing coUimated ^^^Am scans also along side and back surfaces, and by doing 
coUimated single-Compton scattering measurements with higher-energy sources to sample also 
the internal detector volume). 

4.3 Rise time andA/E distributions studies 

A different approach for validating the simulation is given by the comparison of the rise time and 
the A/E distributions as a function of the event energy. Unlike the -^"^^Am measurements which 
create interactions restricted to a region close to the detector surface, 2.6 MeV gamma-rays from 
a ^^^Th source interact in the whole detector volume and can deposit energy at several sites before 
escaping or being absorbed. Under these circumstances, the distributions of the rise times and of 
the pai^ameter A/^ are sensitive to the electric and weighting potential in the whole detector and 
thus are useful for testing the accuracy of the simulation. 

For this study a long measurement was performed with a ^^^Th source^. The detector was 
located in the LNGS underground laboratory and surrounded by a lead shielding in order to reduce 
the background due to cosmic rays and natural radioactivity. 

Figure |T^ shows the distribution of the 10% — 90% rise time as a function of the energy in the 
range including the Double Escape Peak (DEP) at 1592.5 keV and the Full Energy absorption Peak 
(FEP) at 2614.5 keV of ^"^^Tl. Both the experimental and simulated data show a similar- structure 
of the distribution, with a high density band at rise time values around ~ 270 ns. The region below 

^In the simulation only the ^OSjj isotope was considered resulting in some missing gamma-ray interactions at en- 
ergies below '-^1.8 MeV as compared to the experimentally measured spectrum. Moreover, the simulation does not 
include the SEP Doppler broadening of the positron annihilation energy which is determined by the electron momentum 
distribution in germanium. 
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Figure 12. Experimental and simulated distributions of the 10% — 90% rise time values as a function of 
energy from the ^^**Th measurements. In the simulation only ^os^j jg considered. The density plots are 
normalized according to the number of events in the 2.6 MeV FEP. 

the band contains the fast pulses resulting from events with type II charge collection trajectories 
described in Section generated by interactions close to the p+ electrode. The region above the 
band is composed mainly of MSE events which have typically slower rise times than SSE (this 
is visible in Figure Consistently, all the full absorption peaks (consisting typically of multiple 
Compton scatterings followed by photo-absorption) show an important tail in this region, while 
the DEP (containing the typically single-site e^ and e+ absorptions after pair-production) has a 
very weak tail. According to the simulation, events occurring in a small volume close to the outer 
radius of the groove (type III charge collection trajectories in Section are also expected to have 
increased rise times. It is however evident from the plot that the simulation creates events with rise 
times only up to 600 ns while the experimental data show a significant number of events with rise 
times > 600 ns (~ 1.5% of the total). 

The simulated and experimental data can be compared in more detail in Figure |T3| which shows 
distributions of the 10% — 90% rise time in narrow energy regions around the DEP, the Compton 
continuum at 2 MeV and the 2.6 MeV FEP. For the DEP and FEP distributions the contribution 
from their Compton continuum backgrounds is subtracted. While the shape of the distribution 
below ~ 250 ns (corresponding to the fast signals from the type II events) is reproduced fairly 
well, the main peak and the region above it show some differences between the simulation and 
experimental histograms. As we already noted in the previous subsection, our simulation is not 
perfect in some aspects. The charge carrier mobility model is likely responsible for the shift of the 
rise time histogram maximum, while the enhanced drift velocity anisotropy in the simulation could 
explain the greater width of the peak. Minor discrepancies in the geometry description of the p+ 
electrode and the groove, as well as the finite grid width of our pulse library could account for the 
small ~ 20 ns shift between the start of the experimental and simulated histogram, as well as for a 
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Figure 13. Experimental and simulated distributions of the rise time for the region around the DEP, the 
Compton continuum at 2 MeV (CC) and the 2.6 MeV FEP. The error bars take into account only the statistical 
error The DEP and the FEP distribution are corrected for the background. The integrals of the histograms 
are normalized to one. 

part of the excess long rise time events in the experiment. In addition, some pile-up events show 
up with very long rise times up to a few /is. 

However, most of the signals with rise time above ~ 500 ns (most clearly visible in the Comp- 
ton continuum histogram in Figure |T31) are expected to come from the region close to the n+ surface. 
The Li-diffused n-side of the p-n junction can not be fully depleted of conduction band electrons 
and forms the well known dead layer covering the outer surfaces of p-type Ge detectors. The 
appearance of pulses with long rise times in n+ surface interactions was noticed before with con- 
version electron measurements in Ref. lETH and can be clearly observed by irradiating the detector 
with ^'^^ Am gamma rays, as reported in Ref. EBD- It can be assumed that the electric field near the 
n-i- surface has insufficient strength to effectively move the charge carriers. In turn, the drift time is 
extended and part of the charge can be lost by trapping. This effect can explain the observed slow 
pulses in our experimental data. It is also consistent with the fact that the DEP and FEP peak re- 
gions in Figure |T3| show a much smaller excess of slow pulses than the Compton continuum region. 
This is due to the charge carrier losses moving the event in the energy spectrum away from the peak 
regions into the lower energy continuum. As the peak histograms have the continuum background 
subtracted, the slow pulse contribution is removed as well. The content of these slow pulses in the 
Compton continuum region was estimated to be < 4%. The study of these n-i- layer effects can be 
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Figure 14. Experimental and simulated distributions of the AjE parameter as a function of energy for 
^^^Th source measurement. In the simulation only the ^os-pj jg considered. The density plots are normalized 
according to the number of events in the 2.6 MeV FEP. 

in future augmented by our pulse shape simulation if the spatial resolution of the field calculation 
is improved to better cover the 0.7 mm thick region. 

Figure shows the distribution of the parameter A/S as a function of energy. The plots show 
a well defined horizontal band of increased event density, composed of SSE. Below the band the 
region of MSB extends. Above the band one can see a region of events with amplified cuiTcnt signal 
occuiTing close to the p+ electrode as already discussed in Section P3| . This interpretation of the 
A.jE distribution was previously validated by coincident Compton scattering measurements (pro- 
viding clean SSE samples at several energies) and by collimated ^^^Th measurements (providing 
DEP events restricted to the detector volume close and far from the p+ electrode) [0, IBD. Another 
visible feature in the plot is a diagonal band between the DEP and the Single Escape Peak (SEP), 
composed of pair-production events (as explained in Ref. [0, ED), and two weak narrow bands in 
the MSE region, which are caused by cascade summing events'^. There is no significant difference 
between the plots from the experimental and simulated data apart from the ^^^Bi lines present in 
the measurement and not in the simulation. All the features of the A/f" distribution are reproduced 
by the simulation. 

Figure [T3| shows the AjE distribution for energy slices around the DEP, the Compton contin- 
uum at 2 MeV and the SEP at 2103.5 keV. As in Figure |T3| the histograms from the DEP and 
SEP are corrected for the Compton continuum background contribution. The simulation is in good 
agreement with the experimental data in all three energy regions. However, in the Compton con- 
tinuum region one can see a small deviation at low AjE values. This excess can be associated 

'^The two visible additional bands result from the summation of the SSE in the Compton continuum with either the 
511 keV or 583 keV full energy gamma absorption from the ^"^Tl cascade. This results in the energy of the event being 
shifted but the A value staying the same since the SSE from the Compton scattering still dominates the current signal. 
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Figure 15. Experimental and simulated distribution of the AjE parameter in the region around the DEP, 
the Compton continuum at 2 MeV (CC) and the SEP. The error bars take into account only the statistical 
error. The DEP and the SEP distributions are corrected for the background. The integrals of the histograms 
are normalized to one. The top plots show the experimental distribution before and after a cut on the rise 
time to remove slow pulses not included in the simulation. The regions defined by the gray vertical lines are 
enlarged in the bottom plots. 
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to the presence of the near n+ electrode slow pulses in experimental data. It can be assumed that 
by stretching the charge carrier cluster, the peak in the current signal gets wider and its height is 
reduced. By applying a cut on the rise time at 600 ns (black histogram in Figure |T5|) the agreement 
in the tails improves significantly. 

The plots at the bottom of Figure [T3| show the AjE histograms of Compton continuum, the 
SEP and the DEP focused on the region near the SSE-band peak (near A/S ~ 385). It is apparent 
mainly from the DEP plot that the simulated SSE peak is a little wider than in the experimental 
data and the shape is slightly different. This discrepancy is likely the result of the same simulation 
inaccuracies as discussed earlier. The Compton and SEP histograms show that for MSE-dominated 
distributions the overall agreement is significantly better. This demonstrates that the partitioning 
of the energy depositions into individual chai^ge clusters (performed by the Monte Carlo block) is 
simulated accurately. 

5. Background rejection and signal identification studies 

A direct application of the BEGe detector simulation was the study of the Pulse Shape Discrimi- 
nation (PSD) efficiency in separating SSE from MSE introduced in the Section First we com- 
pared the performances of the PSD applied to experimental and simulated sets of measurements 
with ^^^Th and ^°Co sources placed outside the detector housing. We then used the simulation 
to study the survival probabilities of decays which occur inside the Ge crystal. The focus of this 
investigation concerned the Ov/3j3 decay signal of ^^Ge and the background signals created by the 
decays of the radio-isotopes ^^Ge and ^"Co. The latter are produced by cosmic ray interactions 
during germanium detector production above ground. These radio-isotopes are one of the most 
critical backgrounds in ^^Ge neutrinoless double beta decay search. However it is challenging to 
study experimentally the decays of these isotopes inside the Ge crystal. The simulation tools intro- 
duced in this paper represent a powerful method to investigate quantitatively the PSD performance 
on these backgrounds. 

5.1 Comparison of pulse shape discrimination of simulated and experimental data 

The comparison was performed using two high-statistic data sets generated by a ^^^Th source 
^208rpj jj^g simulation) placed 2 cm away from the Al end-cap and a ^"Co source placed at the 
end-cap (on the detector symmetry axis). The recorded signals (10 ns time step) ai^e smoothed 
by triple 50 ns moving window averaging. The cun^ent signal is then reconstructed by 10 ns dif- 
ferentiation of the smoothed charge pulse. Finally, interpolation into 1 ns time steps is performed 
before measuring the current pulse amplitude A. This procedure was found to provide the optimal 
resolution for A. 

The SSE-rich Compton continuum of the 2.6 MeV ^os-pj gamma-line at 2.6 MeV is used to 
calibrate the slight energy dependence of the SSE band in the A /E distribution. This procedure 
was described in detail and validated in Refs. [0, |B1]- The multi-site events, which have lower A/i" 
values than single-site events, are removed by performing the PSD cut at a constant offset from the 
center of the SSE band (visible as a high-density band at A/E' ~ 385 in Figure p^ . 

The efficiency of the cut is tuned by the definition of the offset. In Refs. []7|, ^, the cut offset 
was defined according to the spread of the SSE band in the DEP region and the resulting DEP 
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Table 1. Percentage of survived events from a Th and a Co measurement after the PSD cut in different 
regions of interest (ROI) for the three cut definitions. The resuhs for DEP, SEP, FEP, SP and the two ^"Co 
lines at 1173 keV and 1332 keV (FEPl and FEP2 in the table) are calculated using the net peak areas. The 
CC region corresponds to an 80 keV wide section of Compton continuum centered at 2039 keV (^^Ge Qpi^). 
The uncertainties given in parenthesis for the least significant digits include statistical as well as systematic 
uncertainties estimated by cut parameter variation. The ^os-pj pgp ^j^^ jj^g ^Oq^ FEPl lie outside of the 
range covered by the calibration of A/E energy dependence (1.35 MeV to 2.38 MeV) so can be subject to 
additional systematic errors. 
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acceptance was close to 90%. However, using this definition of the PSD cut on the simulated data 
is expected to create a biased result, since the profile of the SSE-band in our simulation has a 
different shape and spread than the one extracted from the experimental data (see Figure |T3| and 
discussion in the previous section). On the other hand, the profile of the multi-site event region is 
dominated by the partitioning of the charge clusters, so the overestimation of anisotropy or other 
limitations in the simulation have a relatively smaller effect on the MSE as compared to the SSE 
distribution (see Figure |T3|). 

To avoid the effect of the SSE-band discrepancy in the simulation results, we developed al- 
ternative definitions of the A/E cut offset. For experimental data we apply cut offset that retains 
a fixed acceptance of the net area of the DEP at 90%, maintaining compatibility with the previ- 
ously published results. The survival fractions for ^^^Th and ^"Co are listed in the "experimental" 
column of Table |I]. The results are very similar to those obtained in previous works - the small 
differences as compared to []7|] are mainly due to data treatment^ and the results in WITH differ due 
to a different implementation of the analysis. To be able to employ the simulation to estimate a 
realistic experimental survival probability both in SSE-dominated and MSE-dominated data, three 
different versions of the PSD cut were applied to the simulated data. The cut A/E offsets were 
defined by fixing survival fractions in a different peak region for each of the three cuts: 1) in the 

^The different size of the detector used in [0] and geometry of its electrode was found to cause no visible difference 
in the results if the same data treatment is applied. Only the result for the Compton continuum region is significantly 
different, due to a larger distance of the source from the detector, resulting in fewer cascade summing MSE. 
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Figure 16. Experimental (top) and simulated (bottom) energy spectra for a Th source before (red line) and 
after (green line) applying the PSD cut. 



DEP region; 2) in the MSE-dominated SEP of 2°**T1; 3) in the highly-MSE dominated Summation 
Peak (SP) of the two ^"Co gamma-lines. The AjE offsets of these cuts is tuned to be equivalent 
exactly to the experimental data cut with 90% acceptance of DEP events, i.e. the SEP-based cut is 
adjusted to 5.45% SEP survival and SP-based cut to 0.08% SP survival (see the DEP SEP and SP 
entries in the "experimental" column of Table |I]). The results are shown in the last three columns of 
Table H It can be seen that the agreement with the experimental data for Compton continuum^ and 
the MSE dominated FEP and SP regions significantly improves when the SEP-based cut is used 
as compared to the DEP-based cut, while the DEP survival is accordingly underestimated. The SP 
cut results suffer from a higher uncertainty due to the very small number of events remaining in 
the SP after cut, nevertheless the agreement with experimental data is still reasonable in the highly 
MSE-dominated regions. We can conclude that the more similar is the event topology between 
the data region used for cut calibration and the studied region, the better is the agreement between 



simulated and experimental results. Figure [16| shows the simulated and experimental spectra 



before and after applying the PSD cut, using the 90% DEP cut definition. 
5.2 Survival fractions of decays internal the Ge crystal 

In the last part of our study we applied the PSD to simulated ^"Co, ^^Ga decays and OvjSjS decays 
of '^Ge, distributed inside the Ge-detector volume. The aim was to evaluate the OVj8j3-signal 
acceptance by our PSD method and to determine the suppression power of the intrinsic cosmogenic 
backgrounds. Figure |T7| shows the spectra of internal OvjSj3, ^"Co and ^^Ga before and after the 



^Part of the discrepancy in the ^OSjj Compton continuum region arises due to the near n+ layer slow pulses, not 
included in the simulation. 
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Figure 17. Simulated spectra of Ov/3)3 decays (top), ^^Ga (middle) and ^''Co (bottom) inside a BEGe 
detector before (red line) and after (green line) the PSD cuts. 

PSD cut using the DEP-based cut definition. Tlie resulting survival fractions in the energy regions 
around the '^Ge Qpp are summarized in Table ^, including the DEP-based, SEP-based and SP- 
based PSD cut definitions. According to the comparison reported in the previous section, we 
assume that the most realistic value for survival probability is provided by the cut calibrated in the 
data region with the most similar event topology, i.e. the DEP-based cut for the Ovj3j8 peak, the 
SEP-based cut and SP-based cut respectively for the Compton continuum region at 2 MeV of ^^Ga 
and ^''Co (highlighted values). In cases when the MSE content of the region used for calibration 
is higher than in the studied events, the calculated survival probability will be underestimated, and 
vice versa. 

The OVj8j8 spectrum shows a peak at the Qpp energy and a tail extending to low energies 
due to the events in which part of the total energy is either lost in the dead layer or escape the 
detector. Since the electrons generated by OvjSjS decay events have a significant probability to emit 
a bremsstrahlung photon, not all the events in the Qpp peak and in the tail are SSE. The MSE 
contamination is similar to that in the DEP of ^'^^Tl. However, as the energy of the Ovj3j8 decay 
event is higher than in DEP events, also the probability of bremsstrahlung and in turn of an MSE 



Table 2. Percentage of survived events of the simulated Ovj3j3-decays and internal sources of background 
after the cut in an 80 keV wide region centered on the ^^Ge Qpp of 2039 keV. The highlighted numbers 
represent the most realistic predictions for experimental survival probabilities assuming experimental PSD 
cut with DEP fixed to 90%. The uncertainties given in parenthesis for the least significant digits include 
statistical as well as systematic uncertainties estimated by cut parameter variation. The ^"Co result has an 
additional systematic uncertainty due to a small difference on the event topology between the SP region and 
the studied region. 
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interaction is liiglier and the PSD cut survival probability is thus slightly lower. 

The decay chain of the cosmogenic ^^Ge consists of an EC decay into ^^Ga, which then de- 
cays with 1.13 hhalf-Ufe into ^^Zn via EC (11%) and JS"*" decay (89%). The positron has an energy 
endpoint of 1.9 MeV, thus to produce events at the ^^Ge Qjjp(2039 keV) its absorption has to be 
accompanied by an energy deposition from its annihilation photons. Such events have a topol- 
ogy resembling that of SEP events (strong electron/positron interaction combined with a weaker 
511 keV gamma-ray interaction) and thus a comparable survival probability. 

The ^'^Co spectrum shows the two characteristic peaks at 1173.2 and 1332.5 keV generated in 
cascade. As the decay occurs inside the detector, the probability of a coincident detection is high, 
resulting in the summation peak at 2505.7 keV. Moreover, the electron associated to the jS~ decay 
of ^"Co to ^"Ni (endpoint energy of 318 keV) is also absorbed leading to the peculiar "triangular" 
tail at the high-energy side of the peaks. The events in the Qpp region are a result of the cascade 
gamma and beta summation and consequently their topologies are highly multi-site, leading to 
suppression by a factor of 100. 

6. Conclusions and outlook 

In this paper we discussed the development of an integrated simulation tool for the description 
of high-purity Ge detectors and its application to derive the survival probabilities of neutrinoless 
double decays of ^^Ge and of internal decays of ^^Ga and ^'^Co in BEGe type detectors. 

The interaction points of the gamma-rays in the crystal are simulated with the Geant4 based 
MaGe framework; the pulse shapes associated to the interaction points are determined based on 
the electric field calculation with an enhanced version of the MGS program. The outputs of the 
two computations are combined together to determine the final pulse shape corresponding to the 
interaction of the gamma rays within the detector. This pulse shape is then convoluted with the 
front-end electronics response and experimentally measured noise to produce signals that can be 
directly compared to the experimental pulse shapes. 

The particular case of the BEGe detector chosen for the specific application of the simula- 
tion is related to the future use of this type of detectors in the Phase II of the Gerda neutrino- 
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less double beta decay experiment. Previous measurements performed with gamma-ray sources 
have evidenced an enhanced capability of these detectors to discriminate between single-site and 
multiple-site events from the analysis of the signal pulse shape. The characteristic spread in time 
of the signal current pulse created by the different charge clusters of a multi-site event has its origin 
in the sharply rising weighting potential close to the small read-out electrode. 

The validation of the simulation was performed in two different ways. In one of them the 
calculated pulses and their rise times are compared to the measured ones generated in well defined 
positions by the low-energy gamma rays of an ^'^^ Am source. This method is limited to the region 
near the outer surface of the Ge crystal and strongly affected by the definition of the dead layer 
and the size of the spatial grid in the calculation. The other is an integral method and involves 
the comparisons of the rise time and AjE distributions obtained from simulated and experimen- 
tally recorded interactions of high-energy gamma-rays from a 208^^ source with the crystal. Such 
gamma rays interact in the whole volume of the detector and the result of the comparison is a mea- 
sure of the overall goodness of the simulation. Both methods show a good agreement between the 
measured and simulated data. The slight deviations are mainly due to inaccuracy in the geometry 
definition (especially important in case of the p-i- electrode), the finite grid width used for defining 
the library of pulses and the charge carrier mobility model parametrization (measured by using a 
different crystal at a different temperature). 

The simulation was then used to estimate the capability of the pulse shape discrimination 
method to reduce the gamma-ray background produced by internal Ga and ^°Co sources in the 
energy region of the Ovj3j8 decay. The ^^Ge OvjSjS decay was also simulated in order to determine 
the PSD acceptance of its signal. The results are presently dominated by systematic errors related 
to analysis adjustments developed to overcome the imperfections of the simulation, as the AjE 
parameter is particularly sensitive to such effects. Nevertheless, the outcome represents a reason- 
able prediction of the PSD background discrimination efficiency for the Gerda experiment. The 
OVjSjS decay acceptance was found to be (86 it 3)%, while the internal backgrounds were dras- 
tically reduced. Only (0.9+[j 2)% of ^°Co events and (4.5 ±0.3)% of ^*^Ga events at Q^^remain 
after the analysis cut. The ^^Ga suppression can be combined with the time coincidence rejec- 
tion method proposed in O] which takes advantage of the short (1 hour) half life of ^**Ga and the 
10 keV X-ray or Auger electron emission (~ 86% probability) from its mother isotope ^^Ge. This 
latter method, providing reduction by about a factor of 5, will also benefit from the low-noise 
characteristic of BEGe detectors, maintaining a low-energy threshold well below the interesting 
10 keV X-ray energy. Owing to the unique qualities of BEGe detectors, two of the most chal- 
lenging Gerda backgrounds, originating from internal ^°Co and ^^Ga produced by cosmic ray 
interactions in germanium, can be reduced by a factor of 100 by analysis of the germanium signal 
only. 

Acknowledgments 

The authors would like to thank Dino Bazzacco, Bemhard Schwingenheuer, Stefano Riboldi and 
Enrico Borsato for the stimulating discussions and useful comments on this work. They also ac- 
knowledge Matthias Laubenstein for providing help and supporting the measurements in the low 
background laboratory. 



-25- 



This work was supported in part by the Transregio Sonderforschungsbereich SFB/TR27 
"Neutrinos and Beyond" by the Deutsche Forschungsgemeinschaft and by the Munich Cluster of 
Excellence "Origin and Structure of the Universe". 



References 

[1] E. Fiorini, A. PuUia, G. Bertolini, F. Cappellani and G. Restelli, A search for lepton nonconservation 
in double beta decay with a germanium detector, Phys. Lett. B 25 (1967) 602. 

[2] M. Gunther et al Heidelberg-Moscow experiment with Ge-76: Full setup with five detectors, Phys. 
Rev. D 55 (1997) 54. 

C. E. Aalseth et al. (IGEX Collaboration), Recent results oflGEX ^^Ge double-beta decay 
experiment, Phys. of Atomic Nuclei 63 (2000) 1225. 

[3] H.V. Klapdor-Kleingrothaus, A. Dietz, I.V. Krivosheina and O. Chkvorets, Data acquisition and 
analysis of the 76Ge double beta experiment in Gran Sasso 1990-2003, Nucl. Inst. Meth. A 522 
(2004) 371. 

[4] GERDA Collaboration, I. Abt et al., GERDA: The GERmanium Detector Array for the search of 
neutrinoless /3/3 decay of^^Ge at LNGS, Proposal, http://www.mpi-hd.mpg.de/ge76. 

[5] F. Feruglio, A. Strumia and F. Vissani, Neutrino oscillations and signals in /3 and Ovj3/3 experiments, 
Nucl. Phys. B 637 (2002) 345. First KamLAND results, Nucl. Phys. B 659 (2003) 359. 

[6] Canberra Semiconductor NV, Lammerdries 25, B-2430 Olen, Belgium. 

[7] D. Budjas, M. Bamabe Heider, O. Chkvorets, N. Khanbekov and S. Schonert, Pulse Shape Analysis 
with a Broad-Energy Germanium Detector for the GERDA experiment, 2009 JINST 4 P10007. 

[8] D. Budjas, Germanium Detector Studies in the Framework of the GERDA experiment. Dissertation, 
University of Heidelberg (2009). 

[9] M. Barnabe Heider, D. Budjas, K. Gusev and S. Schonert, Operation and performance of a bare 

broad-energy germanium detector in liquid argon, 2010 JINST 5 P10007. 

[10] M. Agostini et al.. Procurement, production and testing ofBEGe detectors depleted in ''^Ge, to appear 
in Nucl. Phys. B, Proc. Suppl.(Neutrino 2010). 

[11] Majorana Collaboration, C. E. Aalseth et al.. The Majorana neutrinoless double-beta decay 

experiment, Phys. Atom. Nucl. 67 (2004) 2002 [Yad. Fiz. 67 (2004) 2025] [arXiv:hep-ex/0405008]. 

[12] P. S. Barbeau, J. I. Collar and O. Tench, Large-Mass Ultra-Low Noise Germanium Detectors: 
Performance and Applications in Neutrino and Astroparticle Physics, JCAP 0709 (2007) 009 
[arXiv:nucl-ex/0701012]. 

[13] Geant4 Collaboration, S. Agostinelli et al. , GEANT4: A simulation toolkit, Nucl. Instrum. and Meth. 
A 506 (2003) 250. 

[14] J. Allison et al., GEANT4 developments and applications, IEEE Trans. Nucl. Sci. 53 (2006) 270. 

[15] M. Bauer et al, MaGe: a Monte Carlo framework for the Gerda and Majorana double beta decay 
experiments, J. Phys.: Conf. Ser. 39 (2006) 362; 

M. Boswell et al, MaGe - a Geant4-based Monte Carlo Application Framework for Low-background 
Germanium Experiments, arXiv:101 1.3827vl. 

[16] P. Medina, C. Santos and D. ViUaume, A simple method for the characterization ofHPGe detectors, 
IEEE IMTC (2004). 



-26- 



[17] L. Mihailescu, W. Gast, R.M. Lieder, H. Brands and H. Jager, The influence of anisotropic electron 
drift velocity on the signal shapes of closed-end HPGe detector, Nucl. Inst, and Meth. A 447 (2000) 
350. 

[18] B. Brayneel, P. Reiter and G. Pascovici, Characterization of large volume HPGe detectors. Part I: 
Electron and hole mobility parametrisation, Nucl. Inst, and Meth. A 569 (2006) 764. 

[19] Zhong He, Review of the Shockley-Ramo theorem and its application in semiconductor gamma-ray 
detectors, Nucl. Inst, and Meth. A 463 (2001) 250. 

[20] CANBERRA Broad Energy Ge (BEGe) Detector, catalogue accessed at URL: 
|http://www.canberra.com/products/485.asp| . 

[21] M. Agostini, Characterization of a Broad Energy Germanium detector through advanced pulse shape 
analysis techniques for tlie GERDA double-beta decay experiment. Master Thesis, University of 
Padova (2009). 

[22] M. Agostini et al., "Characterization of a broad energy germanium detector and application to 
neutrinoless double beta decay search in ^^Ge", companion paper, submitted to JINST, 
arXiv: 1012.5200. 

[23] See |http://www.caen.it/nuclear/product.php?mod=N1728B| . 

[24] See |http://www.iphc.cm-s.fr/-TUC-.htmil . 

[25] M. Agostini, L. Pandola, P. Zavarise and O. Volynets, GELATIO: a general framework for modular 
digital analysis of HPGe signals, in preparation. 

[26] J. Stein, F. Scheuer, W. Gast and A. Georgiev, X-ray detectors with digitalizer preamplifiers, Nucl. 
Inst, and Meth. B 113 (1996) 141. 

[27] M.G. Strauss and R.N. Larsen, Pulse height defect due to electron interaction in the dead layers of 
Ge(li) y-Ray detectors, Nucl. Inst, and Meth. 56 (1967) 80. 

[28] C.E. Aalseth et al.. Results from a Search for Light-Mass Dark Matter with a P-type Point Contact 
Germanium Detector, arXiv:1002.4703v2. 



-27- 



